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£j ■ ABSTRACT 

c/2 \ We report on partially overlapping XMM-Newton (^260 ks) and Suzaku (~100 ks) observa- 

tions of the iron K band in the nearby, bright Seyfert 1 galaxy Mrk 509. The source shows 
a resolved neutral Fe K line, most probably produced in the outer part of the accretion disc. 
Moreover, the source shows further emission blue-ward of the 6.4 keV line due to ionized 
J> ■ material. This emission is well reproduced by a broad line produced in the accretion disc, 

CNj | while it cannot be easily described by scattering or emission from photo-ionized gas at rest. 

The summed spectrum of all XMM-Newton observations shows the presence of a narrow ab- 
sorption line at 7.3 keV produced by highly ionized outflowing material. A spectral variability 
study of the XMM-Newton data shows an indication for an excess of variability at 6.6-6.7 keV. 
These variations may be produced in the red wing of the broad ionized line or by variation of 
a further absorption structure. The Suzaku data indicate that the neutral Fe Ka line intensity 
OS , is consistent with being constant on long timescales (of a few years) and they also confirm as 

■ most likely the interpretation of the excess blueshifted emission in terms of a broad ionized 

Fe line. The average Suzaku spectrum differs from the XMM-Newton one for the disappear- 
ance of the 7.3 keV absorption line and around 6.7 keV, where the XMM-Newton data alone 

^ 1 suggested variability. 

H 
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1 INTRODUCTION 

Deep investigations of the Fe K band in the brightest AGNs allow 
us to probe the presence of highly ionized emitting/absorbing com- 
ponents from the innermost regions around the central black hole. 
The high-sensitivity X-ray satellites XMM-Newton and Chandra 
have shown that the presence of a narrow core of the lowly ionized 
Fe Ka line is nearly ubiquitous (Yaqoob & Padhmanaban 2004; 
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Guainazzi et al. 2006; Nandra et al. 2007) and that ionized compo- 
nents of the line, generally associated with emission from photo- 
and/or collisionally- ionized distant gas are also common (NGC 
5506, NGC 7213, IC 4329A; Bianchi et al. 2003; Page et al. 2003; 
Reynolds et al. 2004; Ashton et al. 2004; Longinotti et al. 2007; 
see also Bianchi et al. 2002; 2005). The presence of broad (neu- 
tral or ionized) components of Fe K lines can only be tested via 
relatively long exposures of the brightest sources (e.g., Guainazzi 
et al. 2006; Nandra et al. 2007). Moreover, the observational evi- 
dence for broad lines and their interpretation in terms of relativistic 
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effects may be questioned when an important absorbing ionized 
component is present. Spectral variability studies help in disentan- 
gling the different, often degenerate, spectral components (Ponti 
et al. 2004; Iwasawa et al. 2004; Ponti et al. 2006; Tombesi et al. 
2007; Petrucci et al. 2007; DeMarco et al., in prep.). 

Mrk 509 (z=0.034397) is the brightest Seyfert 1 of the hard 
(2-100 keV) X-ray sky (Malizia et al. 1999; Revnivtsev et al. 2004; 
Sazonov et al. 2007) that is not strongly affected by a warm ab- 
sorber component (Pounds et al. 2001; Yaqoob et al. 2003). The 
HETG Chandra observation confirms the presence of a narrow 
component of the Fe K line with an equivalent width (EW) of 50 eV 
(Yaqoob et al. 2004). The presence of a second ionized component 
of the Fe K line at 6.7-6.9 keV has been claimed by Pounds et al. 
(2001) who fitted it using a relativistic profile, but Page et al. (2003) 
showed that the same spectral feature was consistent also with a 
simple Compton reflection component from distant material. The 
broad-band BeppoSAX spectrum and, in particular, the soft excess, 
have been fitted by De Rosa et al. (2004) with a reflection com- 
ponent from a ionized disc in addition to a neutral reflection com- 
ponent. Finally, Dadina et al. (2005) found evidence of absorption 
due to transient, relativistically red-blue- shifted ionized matter. 

Here we present the spectral and variability analysis of the 
complex Fe K band of Mrk 509, using the whole set of XMM- 
Newton and Suzaku observations. The paper is organized as fol- 
lows. Section 2 describes the observations and the data reduction. 
In Sect. 3 the spectral analysis of the EPIC-pn data of the Fe K 
band (using phenomenological models) is presented. In particular 
in Sect. 3.3, to check for the presence of an absorption line, the 
EPIC-MOS data have also been considered. In Sect. 3.4 the spec- 
tral variability analysis, within the XMM-Newton observations, is 
presented. Sect. 4 describes the spectral analysis of the Fe K band 
of the Suzaku summed (XIS0+XIS3) data and the detailed com- 
parison with the spectrum accumulated during the XMM-Newton 
observations. In Sect. 4. 1 the HXD-pin data are introduced in or- 
der to estimate the amount of reflection continuum present in the 
source spectrum. Finally, a more physically self-consistent fit of 
the EPIC spectra of all the EPIC instruments (EPIC-pn plus the two 
EPIC-MOS) is investigated in Sect. 5. The results of our analysis 
are discussed in Sect. 6, followed by conclusions in Sect. 7. 



2 OBSERVATIONS AND DATA REDUCTION 

Mrk 509 was observed 5 times by XMM-Newton on 2000-10-25, 
2001-04-20, 2005-10-16, 2005-10-20 and 2006-04-25. All ob- 
servations were performed with the EPIC-pn CCD camera oper- 
ating in small window observing mode and with the thin filter ap- 
plied. The total pn observation time is of about 260 ks. Since the 
live-time of the pn CCD in small window mode is 7 1 per cent, the 
net exposure of the summed spectrum is of about 180 ks. The anal- 
ysis has been made with the SAS software (version 7.1.0), starting 
from the ODF files. Single and double events are selected for the 
pn data, while only single events are used for the MOS camera be- 
cause of a slight pile-up effect. For the pn data we checked that 
the results obtained using only single events (that allow a superior 
energy resolution) are consistent with those from the MOS, finding 
good agreement. The source and background photons are extracted 
from a region of 40 arcsec within the same CCD of the source both 
for the pn and MOS data. Response matrices were generated using 
the SAS tasks RMFGEN and ARFGEN. 

Suzaku observed Mrk 509 four times on 2006-04-25, 2006- 
10-14, 2006-11-15 and 2006-11-27. The last XMM-Newton and 



the first Suzaku observations overlap over a period of ~ 25 ks. 
Event files from version 2.0.6.13 of the Suzaku pipeline processing 
were used and spectra were extracted using XSELECT. Response 
matrices and ancillary response files were generated for each XIS 
using XISRMFGEN and XISSIM ARFGEN version 2007-05- 
14. The XIS 1 camera data are not considered here because of the 
relatively low effective area in the Fe K energy interval, while the 
XIS2 is unavailable for observations performed after November 
2006. We used the data obtained during the overlapping interval 
to check whether the EPIC pn and MOS data on one hand and the 
Suzaku XIS0 and XIS3 data on the other hand are consistent within 
the inter-calibration uncertainties. We found an overall good agree- 
ment between the data from the two satellites, the parameters re- 
lated to the main iron emission features and the power-law contin- 
uum being the same within the errors (except for the XIS2 camera 
above 8 keV). The total XIS observation time is about 108 ks. The 
source and background photons are extracted from a region of 4.3 
arcmin within the same CCD of the source. For the HXD/PIN, in- 
strumental background spectra and response matrices provided by 
the HXD instrument team have been used. An additional compo- 
nent accounting for the CXB has been included in the spectral fits 
of the PIN data. 

All spectral fits were performed using the Xspec soft- 
ware (version 12.3.0) and include neutral Galactic absorption 
(4.2x 10 20 cm -2 ; Dickey & Lockman 1990), the energies are rest 
frame if not specified otherwise, and the errors are reported at the 
90 per cent confidence level for one interesting parameter (Avni 
1976). The sum of the spectra has been performed with the M ATH- 
PHA, ADDRMF and ADDARF tools within the HE AS OFT 
package (version 6.1). 



3 FE K BAND EMISSION OF MRK 509: THE 
XMM-Newton DATA 

The primary goal of this investigation is the study of the Fe K line 
band; therefore, in order to avoid the effects of the warm absorber 
(although not strong; Yaqoob et al. 2003; Smith et al. 2007) and of 
the soft excess, we concentrate on the analysis of the data in the 
3.5-10 keV band only. A detailed study of the warm absorber and 
its variations will be performed by Detmers et al. (in prep), we can 
nevertheless anticipate that the warm absorber has negligible effect 
in the Fe K energy band and thus on the results presented here. 
Figure Q] shows the source light curve in the 3.5-10 keV energy 
band obtained from the XMM-Newton pointings. Mrk 509 shows 
variations of the order of ~30 per cent over the different observa- 
tions, while almost no variability is detected within each observa- 
tion. Only during the fourth observation the source shows signifi- 
cant variability, with a mean fractional rms of about 0.04. 

We start the analysis of the XMM-Newton data considering 
the spectra from the EPIC-pn camera only (including the EPIC- 
MOS data only when a check of the significance of a feature is 
required; see Sect. 3.3). We have fitted a simple power law model 
to the 3.5-10 keV data and found that the spectral index steep- 
ens with increasing flux. It goes from 1.54±0.03 to 1.72±0.03 for 
fluxes of 2.5xl0 -11 and 3.3xl0" n erg cm s , respectively 
(3.0xl0 -11 -4.3x10"" erg cm" 2 s" 1 , in the 2-10 keV band). 
We firstly phenomenologically fitted the Fe K complex of each sin- 
gle observation with a series of emission-absorption lines (see also 
§3.3) and checked that the results on the parameters of Fe K com- 
plex obtained in each observation are consistent within the errors 
(not a surprising result in light of the low statistics of the single 
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Figure 1. 3.5-10 keV EPIC-pn light curves of the XMM-Newton observa- 
tions. The abscissa shows the observation time in seconds. The time be- 
tween the different observations is arbitrary. The black, red, green, blue 
and light blue show the light curves during the 2000-10-25, 2001-04-20, 
2005-10-16, 2005-10-20 and 2006-04-25 observations, respectively. 



spectra and weakness of the ionized features; see § 3.4). Hence, we 
concluded that the continuum variations do not strongly affect the 
observed shape of the narrow-band emission/absorption structures 
in the Fe K band. Thus, in order to improve the signal-to-noise ra- 
tio and thus to detail the fine structures of the Fe K band, the spectra 
of all the XMM-Newton observations have been summed (see t]3.4l 
for the study of the source spectral variability). The summed mean 
EPIC-pn spectrum has been grouped in order to have at least 1000 
counts in each data bin. Moreover, this binning criterion ensures to 
have at least 30 data-points per keV in the 4-7 keV band, where 
the Fe Ka complex is expected to contribute. This guarantees a 
good sampling of the energy resolution of the instrument and the 
possibility of fully exploiting the spectral potentials of the EPIC in- 
struments. Fig. [2] shows the ratio between the data and the best-fit 
power law. The energy band used during the fit has been restricted 
to 3.5-5 and 8-10 keV, in order to avoid the Fe K band, hence mea- 
suring the underlying continuum. The resulting best-fit power-law 
continuum has a photon index of 1.63±0.01 and very well repro- 
duces the source emission (X 2 =170.0 for 163 de erees of freedom, 
dof) outside the Fe K band. The inclusion of the Fe K band shows 
that other components are necessary to reproduce it (x 2 =753.9 for 
307 dof). The bad statistical result is explained by the presence of 
clear spectral complexity in the 6-7 keV band. 



3.1 The 6.4 keV emission line 

Panel a of Fig.|2]shows the clear evidence for a prominent emission 
line, consistent with a neutral Fe Ka line at 6.4 keV. We there- 
fore added a Gaussian emission line to the model, obtaining a very 
significant improvement of the fit (Ax 2 =392.1 for the addition of 
3 dof). The best-fit energy of the line is 6.42±0.02 keV, consis- 
tent with emission from neutral or slightly ionised material. The 
line has an equivalent width of 69±8 eV and is clearly resolved 
(cr=0.12±0.02 keV), as shown by the contour plot in the left panel 
of Fig. [3] The residuals in panel b of Fig.|2]show no excess redward 
of this emission line, which could have been indicative of emission 
from relativistically redshifted neutral material. 
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Figure 2. (Upper panel) Observed-frame 3.5-10 keV summed XMM- 
Newton EPIC-pn spectrum fitted in the 3.5-5 and 8-10 keV band with a 
power law. (Panel a) Data/model ratio. This ratio shows a clear evidence 
for a neutral Fe K emission line and further emission from ionized Fe, as 
well as other complexities around 7 keV. (Panel b) Data/model ratio when 
two resolved emission lines (for the Fe Ka and K/3) are included in the 
spectral fitting. Strong residuals are still present, indicative of ionized Fe K 
emission, while no residual emission redward of the neutral Fe K line ap- 
pears. (Panel c) Data/model ratio when a narrow emission line is included 
in the model to reproduce the ionized emission. (Panel d) Same as panel c, 
but with a single broad emission line instead of a narrow line. In both cases 
(panel c and d), an absorption feature is present around 7 keV. (Panel e) 
Data/model ratio when an absorption component (modeled using XSTAR) 
and a relativistic ionized line are added to the power law and the emission 
from neutral Fe K. 
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Figure 3. (Left panel) Contour plot of the sigma vs. intensity of the neutral 
Fe K line. (Right panel) Contour plot of the energy vs. intensity of the nar- 
row line used to fit the ionized Fe K emission. The narrow line energy is not 
consistent with emission from Fe XXV (neither with the forbidden at 6.64 
keV, nor with the resonant at 6.7 keV), Fe XXVI or Fe K/3, whose energy 
is indicated by the vertical dotted lines (from left to right). 



3.2 The ionized Fe K emission line 

An excess is, however, present in the range 6.5-7 keV (Fig.[2] panel 
a). If modeled with a Fe K/3 component with the expected energy 
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(fixed at 7.06 keV) and forced to have an intensity of 0.15 of the 
Ka (Palmeri et al. 2003a,b; Basko 1978; Molendi et al. 2003) and 
a width equal to the Fe Ka line (i.e. assuming that the Ka and K/3 
line originate from one and the same material), the fit improves sig- 
nificantly (Ax 2 =20.3). Nonetheless, significant residuals are still 
present in the 6.5-6.9 keV band (panel b of Fig(2]l. If this further 
excess is modelled with a narrow Gaussian line (Ax 2 =25 for 2 
additional dof), the feature (EW=12±4 eV) is found to peak at 
E=6.86±0.04 keV (see panel c of Fig. 2 and right panel of Fig. 
0. Thus, the energy centroid is not consistent with the line being 
produced by either Fe XXV or Fe XXVI (right panel of Fig. [3} in 
a scattering medium distant from the X-rays source (Bianchi et al. 
2002; 2004). The higher energy transition of the Fe XXV complex 
is the "resonant line" expected at 6.7 keV (see e.g. Bianchi et al. 
2005). Thus, to save this interpretation, it is required that the photo- 
ionized gas has a significant blueshift (~5700 km/s, if the line is as- 
sociated to Fe XXV) or redshift (~4500 km/s, for Fe XXVI). Then, 
instead of fitting the ionized excess with a single line, we fitted it 
with two narrow lines forcing their energies to be 6.7 and 6.96 keV. 
The fit clearly worsens (x 2 =326.7 for 302 dof, corresponding to a 



Ax 2 



-10.1 for the same dof). However, if the gas is allowed to 



be outflowing, the fit improves (Ax 2 =4.3 for the addition of 1 new 
parameter; x 2 =312.3 for 301 dof; the EW are 8.9 and 12.4 eV for 
the Fe XXV and Fe XXVI lines, respectively) as respect to the sin- 
gle narrow emission line and it results to have a common velocity 

of 35001 "oo km/s- 

Alternatively, the excess could be produced by a single broad 
line coming from matter quite close to the source of high-energy 
photons (in this case the Fe K emission is composed by Fe Ka+/3 
plus another Fe K line). Leaving the width of the line free to vary, 
the fit improves, with x 2 =311.1 (panel d Fig. 2) and A^ 2 of 5.5, 
with respect to the single narrow ionized emission line fit, and A^ 2 
of 1.3 for the same dof with respect to the best-fit model with two 
narrow ionized lines. The resulting broad ionized Fe K line has 
EW=23±9 eV and cr=0. 14±^;J| keV. The best-fit energy of the line 
does not change significantly (E=6.86to' ( j , 5 keV); however, in this 
case the emission is consistent (at the 99 per cent confidence level) 
with either Fe XXV or Fe XXVI. Although the statistical improve- 
ment is not highly significant, in the following we will consider 
that the ~6.8-6.9 keV excess is indeed associated with a resolved 
emission line. 



3.3 Ionized absorption? 

The XMM-Newton data also display a narrow absorption feature at 
E~7 keV (observed frame; see Fig. 2, panel d). Since this feature is 
very close to the broad excess we just discussed, its significance and 
intensity are degenerate with the broad emission-line parameters. 
In order to gain some insight, we then fixed the broad emission- 
line parameters at the best-fit ones obtained before the addition of 
a narrow (a fixed at 1 eV) Gaussian absorption line component. In 
this case, the line is significant at the ~99 per cent confidence level 
(dashed contours of Fig.|4j Ax 2 =15.5 for 2 additional parameters; 
see also panel e of Fig. 2). Once the MOS dat£] are added, the 
significance of this feature increases to 99.9 per cent (solid contours 
of Fig.[4ji, in both cases, of a broad and of a narrow ionized emission 
line. The best fit energy and EW of the line are E=7.28 ±° ™ keV 




observed energy (keV) 

Figure 4. Superposition of the pn (green), MOS1 (black) and MOS2 (red) 
summed spectra of all the XMM-Newton observations. The data are fitted, 
in the 3.5-5 and 8-10 keV bands with a power law, absorbed by Galactic 
material. The same structures are present in the three spectra. In particular, a 
narrow drop of emission is present in all the instruments at the same energy 
(see vertical dotted line). (Inset panel) Confidence contour plot of the inten- 
sity vs. energy of the narrow unresolved ionized absorption when using the 
pn data alone (dashed contours) and including the MOS data as well (solid 
contours). The lines indicate the 68.3 (black), 90 (red), 99 (green) and 99.9 
(blue) per cent confidence levels. 



andEW=-14.9^; 2 eV, E=7.33±^; ^ keV andEW=-13.ll|| eV, 
in the pn alone and in the pn+MOS, respectively. 



3.4 Time resolved spectral variability and total rms spectrum 

One of the goals of the present analysis is to search for time- 
variation of the emission/absorption features of the Fe K complex. 
To measure possible variations in the Fe K band, the mean EPIC- 
pn spectra of each of the 5 XMM-Newton observations have been 
studied. The spectra are fitted with the same model composed by 
a power law plus three emission lines for the Fe Ka, K/3 (with the 
width fixed at the best-fit value, a=0. 12 keV) and the broad ionized 
Fe K line. The low statistics of the spectra of the single observations 
prevents us from the detection of significant spectral variability of 
the weak ionised emission/absorption lines. The neutral Fe K line 
is better constrained and we find that its EW is anti-correlated with 
the level of the continuum, as expected for a constant line. 

A different, more sensitive, way to detect an excess of spectral 
variability is the total rms function. The upper panel of Fig.[5]dis- 
plays the shape of the summed spectrum in the Fe K line band. The 
lower panel shows the total rms spectrum (Revnivtsev et al. 1999; 
Papadakis et al. 2005) calculated with time bins of -~4.5 ks. The 
total rms is defined by the formula: 



RMS(E) 



AE * arf(E) 



(1) 



where S 2 is the source variance in a given energy interval AE; 
< <r 2 rl , > is the scatter introduced by the Poissonian noise and art 
is the telescope effective area convolved with the response matrijfl 



The shapes of the emission/absorption lines in the MOS instruments 
appear slightly narrower, although consistent with the values obtained with 
the pn instrument. 



2 The total rms spectrum provides the intrinsic source spectrum of the 
variable component. Nevertheless, we measure the variance as observed 
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This function shows the spectrum of the varying component only, 
in which any constant component is removed and has been com- 
puted by using the different XMM-Newton observations as if they 
were contiguous. The total rms spectrum may be reproduced by a 
power law with a spectral index of 2. 13 (x 2 =46.4 for 43 dof). Thus, 
the variable component is steeper than the observed power law 
in the mean spectrum, in agreement with the mentioned observed 
steepening of the photon index (F) with flux. The T-fiux correla- 
tion is commonly observed in Seyfert galaxies and has been inter- 
preted as being due to the flux-correlated variations of the power- 
law slope produced in a corona above an accretion disc and related 
to the changes in the input soft seed photons (e.g. Haardt, Maraschi 
& Ghisellini 1997; Maraschi & Haardt 1997; Poutanen & Fabian 
1999; Zdziarski et al. 2003). These models predict the presence of 
a pivot point, that would correspond to a minimum in the total rms 
spectrum. The observation of a perfect power law shape (see Fig. 
5) indicates that the pivot point (if present) has to be outside the 
3-10 keV energy band. On the other hand, the slope-flux behaviour 
can be explained in terms of a two-component model (McHardy, 
Papadakis & Uttley 1998; Shih, Iwasawa & Fabian 2002) in which 
a constant-slope power law varies in normalization only, while a 
harder component remains approximately constant, hardening the 
spectral slope at low flux levels only, when it becomes prominent 
in the hard band. In this scenario the spectral index of the variable 
component is equal to the one of the total rms spectrum, that is 
F=2.13. 

Moreover we note that at the energy of the neutral and ionized 
Fe K line components, no excess of variability is present, in agree- 
ment with these components being constant, while an indication for 
an excess of variability is present around 6.7 keV. In order to com- 
pute the significance of this variability feature, a narrow Gaussian 
line has been added to the modelling of the total rms spectrum. The 
best-fit energy of the additional line is 6.69 keV, with a a fixed 
at the instrumental energy resolution, while the resulting A\ 2 is 
8.9 for the addition of 2 parameters (that corresponds to an F-test 
significance of 98.8 per cent). Introducing the line, the continuum 
spectral index steepens to F ~ 2.18. The dashed line in Fig. 5 high- 
lights the centroid energy of the neutral Fe Ka line, while the dotted 
line (at ~6.7 keV, rest frame) is placed at the maximum of the vari- 
ability excess. This energy corresponds to a drop of emission in the 
real spectrum, as we shall discuss in more detail in Section 5. 



4 THE Suzaku VIEW OF THE FE K BAND EMISSION 

As mentioned in §2, the source was also observed with Suzaku. 
The first 25 ks Suzaku observation is simultaneous with the last 
XMM-Newton pointing. The source spectra of all the instruments 
are in very good agreement, during the simultaneous observation. 
The spectrum is also consistent with the presence of the emission 
and absorption lines, as observed in the mean XMM-Newton spec- 
trum, nevertheless, due to the low statistics of the 25 ks spectrum 
and the weakness of the ionized features, it is not possible to per- 
form a detailed comparison. Only the presence of the strong Fe Ka 
line can be investigated, the ionized emission and absorption lines 
are not constrained in the 25 ks Suzaku exposure. 

through the instrument. Thus, the sharp features in the source spectrum, 
as well as the effects of the features on the effective area, are broadened by 
the instrumental spectral resolution. For this reason, to obtain the total rms 
spectrum, we take into account the convolution of the effective area with 
the spectral response. 
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Figure 5. Lower panel: Total rms variability spectrum of the XMM-Newton 
observations. The data (blue crosses) show the spectrum of the variable 
component. The best-fit model is a power law with spectral index T=2.18 
(red line) plus a Gaussian emission line (improving the fit by A\ 2 of 8.9 for 
the addition of 2 parameters). The dashed line highlights the centroid energy 
of the neutral Fe Ka line, while the dotted line is placed at the maximum of 
the variability excess, modeled with the Gaussian emission line. The excess 
variability energy corresponds to a drop of emission of the real spectrum. 



Also during the 4 Suzaku pointings, Mrk 509 has shown little 
variability, with flux changes lower than 10-15 per cent, hamper- 
ing any spectral variability study. Fig. [6] shows the XMM-Newton 
(black) and Suzaku XIS0+XIS3 (red) summed mean spectra. The 
data were fitted, in the 3.5-5 and 7.5-10 keV bands, with a simple 
power law and Galactic absorption: the ratio of the data to the best 
fit model is shown in Fig. 6. The source emission varied between 
the XMM-Newton and the Suzaku observations. The best-fit spec- 
tral index and the 3.5-10 keV band fluxes are: r=1.63±0.01 and 
r=1.71±0.02 and 2.63xlO" n and 3.11xl0~ n ergs cm" 2 s _1 , 
during the XMM-Newton and Suzaku observations, respectively. 
The neutral and ionized Fe K emission lines appear constant, while 
some differences are present at 6.7 keV, the same energy where 
the XMM-Newton data were suggesting an increase of variability. 
Other more subtle differences appears at ~ 7 keV, where the ab- 
sorption line imprints its presence in the XMM-Newton data only. 

The Suzaku spectrum of Mrk 509 shows, in good agreement 
with the XMM-Newton one, a resolved neutral Fe K line smoothly 
joining with a higher energy excess, most likely due to ionized iron 
emission (see Fig. |6). Given that no absorption lines around 6.7 
keV or 7.3 keV are present in the Suzaku data, the spectrum may 
be useful to infer the properties of the emission lines more clearly. 

The XIS0+XIS3 Suzaku summed spectrum has been fitted in 
the 3.5-10 keV band with a power law plus two resolved Gaussian 
emission lines to reproduce the emission from Fe Kq+/3. The pa- 
rameters of the Fe Ka line are free to vary, while the Fe K/3 ones 
are constrained as in §3.2. This fit leaves large residuals (x 2 =1379.8 
for 1337 dof) in the Fe K band. In this respect, it is difficult to de- 
scribe the >6.5 keV excess with a single narrow ionized Fe line 
(either due to Fe XXV or Fe XXVI). In fact, although the addition 
of a narrow line is sig nificant (Ax 2 =20.1 for 2 more parameters), 
it leaves residuals in the Fe K band. This remaining excess can 
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3.5-10 keV 



BEST-FIT 



SPECTRA 



A 
B 



Suzaku 

r 

1.72±0.02 
1.72±0.02 



pi norm a 

1.12±0.02 
1.12±0.02 



E Neut. 

keV 
6.42±0.03 
6.42±0.02 



°"Neut. 

keV 
<0.06 
<0.07 



^Neut. b (EW) C 

1.7±0.5 (32) 
2.1±0.5 (40) 



Elon. 

keV 
6.54±0.09 
6.61±0.08 



<?lanJ T in 
keV/r g 

0.40±0.1 
24±10 



Aion. b (EWf 

4.6±1.2 (90) 
3.4±0.8 (79) 



X 2 /dof 

1344/1340 
1346/1340 



Self-consistent model 
XMM-Newton 

T pi norm a 



1.70±0.01 

5 8+ 5 ' 2 

J '°-4.8 



0.92±0.04 

log«) 

5 15+ 1 ' 25 

J - 1 -0.52 



ENeut. 

keV 
6.41±0.01 

z 

-0.0484+2°3 



°"Ncut. 
keV 
0.07±0.01 

X 2 /dof 
894.3/876 



A N cut. 

2.2±0.3 



Incl. 

deg 
47±2 



erg cm s 
,+200 
-7 



ll" 1 



A Refl.Ion. 



0.9 



+3.0 



Table 1. Top panel: Best-fit values of the summed spectra (XIS0+XIS3) of all Suzaku observations fitted in the 3.5-10 keV band. Both model A and B include 
a power law and two Gaussian lines Kct+/3 to fit the 6.4 keV excess. In addition to this baseline model, either another Gaussian component (Model A) or a 
DISKLINE profile (Model B) have been added to reproduce the ionized line, respectively. In Table the best-fit power law spectral index (I 1 ) and normalization 
as well as the Fe Kct energy, width and normalization are reported for model A and B. The energy, width and normalization are reported when a Gaussian 
profile for the ionized Fe K line is considered (Model A), while the best fit energy, inner radius and normalization are presented when a DISKLINE profile 
is fitted (Model B). Standard disc reflectivity index, outer disc radius and disc inclination of a = 3, v ut = 400 r 9 and 30° have been assumed for the 
relativistic profile. Bottom panel: Best fit results of the summed XMM-Newton EPIC-pn and EPIC-MOS data of Mrk 509 (fitted in the 3.5-10 keV band). The 
model (WABS*ZXIPCF*(POW+ZGAUS+ZGAUS+PEXRAV+KDBLUR*(REFLION))) consists of: i) a power law; ii) two Gaussian emission lines for the Fe Ka 
and K/3 emission (this latter has energy is fixed to the expected value, 7.06 keV, intensity and width tied to the Ka values); iii) a neutral reflection continuum 
component (PEXRAV in Xspec) with ii=l (value broadly consistent with the pin constraints and the values previously observed; De Rosa et al. 2004), Solar 
abundance and high energy cut off of the illuminating power law at 100 keV; iv) a ionized disc reflection spectrum (reflion model; Ross & Fabian 2005) with 
the disc inner and outer radii and the emissivity of 6, 400 r 9 and —3, respectively. The best fit disc inclination and ionization and the normalization of the 
disc reflection component are shown; v) an ionized absorption component (ZXIPCF) totally covering the nuclear source. The best fit column density, ionization 
parameter and outflow velocity are reported, a) In units of 10~~ 2 photons keV~ x cm~ 2 s —1 at 1 keV; b) In units of 10~ 5 photons cm~ 2 s —1 ; c) In units ofeV; 
d) In units of 10 22 atoms cm~ 2 . 




Energy (keV) 

Figure 6. XMM-Newton (black) and Suzaku XIS0+XIS3 (red) summed 
mean spectra. The data are fitted, in the 3.5-5 and 7.5-10 keV bands, with 
a simple power law, absorbed by Galactic material, and the ratio of the data 
to the best fit model is shown. The arrows mark absorption features in the 
spectrum. 

be reproduced (Ax 2 =5.9 for 1 more parameter), in a photoionized 
gas scenario, by a blend of two unresolved ionized lines, requir- 
ing three emission lines to fit the Fe K band (FeKa+/3, Fe XXV 
and Fe XXVI). In this case, such as in the analysis of the XMM- 
Newton mean spectrum, a blueshift of this component is suggested 
(v=2600t 2 i )( ) km s _1 ). However, the best-fit model (this scenario 
is strengthened by the lack of narrow peaks) suggests that the ex- 
cess may be in fact associated with a broad ionized Fe line (over 
which the ~ 6.7 keV and ~ 7.3 keV absorption lines are most likely 
superimposed, but during the XMM-Newton observation only). In 
fact considering a broad Fe line instead of the two narrow lines we 



obtain an improvement of Ax 2 =9.7 for the same dof (see TableQ] 
model A). 

Thus, the Suzaku data indicate that the broad excess at 6.5- 
6.6 keV is indeed due to a broad line rather than a blend of narrow 
ionized Fe lines. Since broad lines may arise because of relativis- 
tic effects in the inner regions of the accretion flow, we tested this 
hypothesis by fitting the excess at 6.5-6.6 keV with a DISKLINE 
profile. The statistics of the spectrum is not such to allow us to con- 
strain all the parameters of the ionized DISKLINE model. Thus, the 
disc reflectivity index has been fixed at the standard value (a = —3, 
where the emissivity is proportional to r a ), the outer disc radius 
and inclinations to 400 gravitational radii (r 9 ) and 30°, respectively. 
The broad line is consistent with being produced in the accretion 
disc (Table 1, Model B); however, the emission from the innermost 
part of the disc is not required, the lower limit on the inner disc 
radius being 10-15 r 9 . As clear from Fig. [6] the Suzaku data do not 
require any ionized Fe K absorption structures. 

In order to quantify the differences between the Suzaku and 
XMM-Newton spectra (and, in particular, the reality of the absorp- 
tion structures at 6.7 and 7.3 keV appearing in the XMM-Newton 
spectrum only) we fixed all the parameters of the Suzaku model 
(apart from the intensity and spectral index of the direct power law) 
and fit the XMM-Newton data with that model. This corresponds 
to assuming that the intrinsic line shapes do not vary between the 
two observations. Then, a narrow Gaussian line has been added to 
the XMM-Newton model to estimate the significance of the puta- 
tive absorption structures. The improvement in the spectral fitting 
is evident, as indicated by the Ax 2 =28.3 and 22 in the case of a 
line at E=6.72±0.04 keV and E=7.29±0.04 keV respectively. The 
presence of these spectral features only in the XMM-Newton ob- 
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servations is thus indicative of variability at energies ~6.6-6.7 and 
~7.3 keV. 

4.1 The Suzaku pin data to constrain the reflection fraction 

We add the pin data to measure the amount of reflection contin- 
uum. We note that the pin data provide a good quality spectrum 
up to 50 keV. The model used involves a direct power law plus a 
neutral reflection component (pexrav model in Xspec; Magdziarz 
& Zdziarski 1995) plus the Fe Ka+/3 resolved lines and a broad 
(DISKLINE) component of the line. As for model B we fix some 
of the parameters of the DISKLINE profile (disc inclination=30°, 
r out =400 i g and a=-3). Moreover we assume a high-energy cut off 
of 100 keV and Solar abundance. Thus, by fitting the 3-50 keV 
band data, we obtained a reflection fraction R — OA^g ® and a 
spectral index F=1.76lo 03- The total EW of the emission lines 
above the reflected continuum (about 1.2 keV) is broadly consis- 
tent with the theoretical expectations (Matt et al. 1996) and with 
what observed in Compton thick Seyfert 2 galaxies, where the pri- 
mary continuum is absorbed and only the reflection is observed. 
Nevertheless, also for this source, as already known from previous 
studies (Zdziarski et al. 1999), we observe that the spectral index 
and the reflection fraction are degenerate and strongly depend on 
the energy band considered. In fact, if the 2-10 keV band is con- 
sidered, the reflection fraction increases, resulting to be _R=1.1+q' 2 
and the power-law photon index of r=1.88to Q2- The total EW of 
the Fe emission lines above the reflected continuum are about 750 
eV. Again these values are broadly in agreement with expectations 
(Matt et al. 1996). 



5 A PHYSICALLY SELF-CONSISTENT FIT: POSSIBLE 
ORIGIN OF THE SPECTRAL FEATURES 

The analysis of the XMM-Newton and Suzaku data shows evidence 
for the presence of: i) a resolved, although not very broad, (er ^0.12 
keV) neutral Fe Ka line and associated Fe K/3 emission; ii) an ion- 
ized Fe K emission line inconsistent with emission from a distant 
scattering material at rest and most likely produced in the accretion 
disc; iii) an absorption line at ~7.3 keV, present in the summed 
spectrum of all XMM-Newton observations only; iv) an indication 
for an enhancement of variability - both by considering the XMM- 
Newton data alone and by comparison between the two data sets - at 
~6.7 keV that could be either due to the high variability of the red 
wing of the broad ionized Fe K line, possibly associated with a vari- 
ation of the ionisation of the disc, or to a second ionized absorption 
line. 

These emission/absorption components are partially inter- 
connected to each other given the limited CCD resolution onboard 
XMM-Newton and Suzaku. Thus we re-fit the XMM-Newton (both 
the pn and MOS in the 3.5-10 keV energy band) data with a model 
containing components that better describe the physical processes 
occurring in the AGN. In particular, we consider two Gaussian lines 
for the Fe Ka and K/3 emission plus a neutral reflection compo- 
nent (PEXRAV in XSPEC) with a reflection fraction R = 1 (con- 
sistent with the constraints given by the Suzaku pin data). The Fe 
Ka line has an equivalent width of 1 keV above the reflection 
continuum. Moreover, we fit the broad ionized Fe K line with a 
fully self-consistent relativistic ionised disc reflection component 
(reflion model in Xspec; Ross & Fabian 2005, convolved with a 
LAOR kernel; KDBLUR in Xspec). 

The statistics prevents us from constraining the parameters of 



the relativistic profile. Standard values for the relativistic profile are 
assumed, with the disc inner and outer radii and the emissivity of 6, 
400 r 9 , and —3, respectively. Finally, the ~7.3 keV absorption line 
has been fitted with a photoionised absorption model (zxipcf model 
in Xspec; Miller et al. 2007; Reeves et al. 2008; Model C, Table 1), 
assuming a total covering factor. 

Table [TJ shows the best-fit parameters. Once the presence of 
the reflection continuum is taken into account, the power law slope 
becomes steeper (F=1.70±0.01, AF ~0.07) as compared to the fit 
with a simple power law and emission absorption lines (see §3). 
The best fit energy of the neutral Fe Ka line is E=6.41±0.01 keV, 
consistent with being produced by neutral material, and results to 
be narrower (<r=0.07±0.01 keV) than in the previous fits. The ion- 
ized emission line is fitted with a ionized disc reflection model. 
The only free parameters of such a component are the inclina- 
tion and ionisation parameter of the disc that result to be 47±2° 
and £=ll+;f erg cm s" 1 (Model C, Table 1). The material pro- 
ducing the 7.3 keV absorption feature in the XMM-Newton data 
has to be highly ionized, as also indicated by the absence of a 
strong continuum curvature. In fact, the best ionization parameter 
is log(0=5.15t o 2 5 5 2 and the column density N H = 5.8t 5 4 2 s x 10 22 
cm~ 2 . Nevertheless the observed energy of the absorption feature 
does not correspond to any strong absorption features, thus there 
is evidence for this absorption component to be outflowing with a 
shift v = -0.0484t°;^ 2 c (~ 140001^° km s _1 ). The resulting 
X 2 is 894.3 for 876 dof. 



6 DISCUSSION 

This study clearly shows that long exposures are need to disentan- 
gle the different emitting/absorbing components contributing to the 
shape-variability of the Fe K complex in Seyfert galaxies. Here 
we discuss the origin of both neutral and ionized emission and ab- 
sorption Fe lines in Mrk 509 which allow to have insights in the 
innermost regions of the accretion flow. 



6.1 Neutral/lowly ionized Fe emission line 

Once the broad ionized line is fitted, the width of the Fe Ka line 
lowers to a value of 72±11 eV (see Fig. [3} that corresponds to 
a FWHMfFe Ka)=8000±1300 km s _1 (see Model C, Table 1). 
This value is slightly higher than that measured by Yaqoob & Pad- 
manabhan with a ~50 ks HETG Chandra observation (2820^2800 
km s _1 ). The FWHM of the Fe Ka line is larger than the width 
of the H/3 line (FWHM(H/3)=3430±240 km s" 1 ; Peterson et al. 
2004; Marziani et al. 2003), indicating that the Fe line is produced 
closer to the center than the optical BLR and, of course, than the 
torus postulated in unified models; we note that a wide range of 
FWHM values is observed for the BLR and the Fe K lines in local 
Seyfert galaxies (Nandra 2006). However, the UV and soft X-ray 
spectra of Mrk 509 show evidence for the presence of broad emis- 
sion lines with FWHM of 11000 km s" 1 (Kriss et al. 2000). The 
origin of these UV and soft-X lines is still highly debated, nev- 
ertheless they may indicate that the BLR region is stratified, i.e. 
that these lines are not produced in the optical BLR but in the in- 
ner part of a stratified BLR region (see also Kaastra et al. 2002; 
Costantini et al. 2007), possibly as close as 2000 r 9 from the cen- 
ter (about 0.012 pc, being the mass of the black hole in Mrk 509 
Mbh ~1.43±0.12xl0 8 M q Peterson et al. 2004; Marziani et al. 
2003). Nevertheless, if the line is produced in the innermost part 
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of a stratified BLR, it would require either a higher covering frac- 
tion or a higher column density than generally derived from the 
optical and ultraviolet bands. Simulations by Leahy & Creighton 
(1993) show that about 70 per cent of the sky, as seen by the cen- 
tral source, has to be covered in order to produce the Fe Ka line, 
if the broad line clouds have column densities of about 10 23 cm -2 , 
while the typical values for the BLR clouds covering fractions are 
of the order of 10-25 per cent (Davidson & Netzer 1979; Goad & 
Koratkar 1998). Alternatively, the Fe Ka line may be produced by 
reflection by the outer part of the accretion disc. 



6.2 Ionized Fe emission lines 

The spectrum of Mrk 509 shows emission from ionized iron, con- 
sistent with either Fe XXV or Fe XXVI, implying photoionized gas 
outflowing or inflowing respectively. Alternatively, the ionized Fe 
K emission may be produced by reflection from the inner part of 
the accretion disc. 

In fact, both the XMM-Newton and the Suzaku data are consis- 
tent with the two scenarios, even if a slightly better fit (Ax 2 =5.5 
and 9.7 for XMM-Newton and Suzaku, respectively) is obtained 
in the case of broad line. Moreover in the case of narrow emis- 
sion lines the emitting gas should have a significant outflow (for 
Fe XXV, v-3500 and 2600 km s _1 for XMM-Newton and Suzaku, 
respectively) or inflow (for Fe XXVI, v~4500 km s _1 ) with ve- 
locities higher than what generally observed (Reynolds et al. 2004; 
Longinotti et al. 2007; but see also Bianchi et al. 2008 that detect 
an outflow of v=900±7qo km s _1 ). On the other hand, the high 
radiative efficiency of the source (77=0.12; Woo & Urry 2002) sug- 
gests that the accretion disc is stable down to the innermost re- 
gions around the BH, where the reflection component should be 
shaped by relativistic effects. For these reasons, although an out- 
flowing emitting gas is not excluded, the broad line interpretation 
seems favoured. In fact, the profile of the line is compatible with 
being shaped by relativistic effects, consistent with its origin be- 
ing in the surface of an accretion disc, in vicinity of a black hole. 
Nevertheless, the width of the line is not a compelling evidence. 
The observed broadening of the line can be reproduced also with 
the Comptonization process occurring in the upper layer of the ion- 
ized accretion disc. Moreover, we stress that the main evidences for 
the presence of a broad Fe K line comes from the mean summed 
spectrum. The process of summing spectra, although is a power- 
ful way to extract information, might be dangerous in presence of 
spectral variability and when applied to observations taken many 
years apart. Thus, the final answer on the origin of these ionized 
lines will be obtained with either a higher resolution observation or 
with significantly longer XMM-Newton exposures. 

6.3 Ionized Fe absorption lines 

The XMM-Newton data indicate the presence of a highly ion- 
ized absorption component, the best fit column density being 
N H =5.8t4 gXl0 22 cm' 2 and ionization log(0=5.15tJ g. More- 
over, fitting the absorption with this model, it results that the ab- 
sorber has to be blueshifted by 0.0484+° °i 2 c. The blueshift corre- 
sponds to an outflow velocity of ~14000 km s - . The structure im- 
plies a significant blueshift if the absorber is located in the core of 
Mrk 509 but, considering the systemic velocity of the galaxy, its en- 
ergy is also consistent with a local absorber (McKernan et al. 2004; 
2005; Risaliti et al. 2005; Young et al. 2005; Miniutti et al. 2007; 
but see also Reeves et al. 2008). Nevertheless, the observed vari- 



ability between the XMM-Newton and Suzaku observations points 
towards an origin within Mrk 509. 

An hint of variability is observed around 6.7 keV both in 
the XMM-Newton data and by comparing the XMM-Newton and 
Suzaku spectra. This could be due in principle to variability in 
the red wing of the ionized emission line. However, the total rms 
spectrum shows a peak of variability that is consistent with be- 
ing narrow, thus it may suggest an alternative explanation. Indeed, 
the observed difference between the XMM-Newton and Suzaku Fe 
K line shapes could be due to a further ionized absorption com- 
ponent, present only the XMM-Newton observations, with a col- 
umn density Nh=5.4_ 4 4 x 10 cm and ionization parameter 
log(O=2.04i ( ( ];gf ) . When the structure at 6.7 keV is fitted with such 
a component, an absorption structure appears around 7.3 keV, nev- 
ertheless its equivalent width is not strong enough to reproduce the 
total absorption feature; moreover, it appears at slightly different 
energy, not completely fitting the ~7.3 keV line. Thus, the absorp- 
tion structures at 6.7 and the one at 7.3 keV may be connected 
and they may be indicative of another absorption screen. If this 
further lower ionization absorption component is present, different 
absorption feature would be expected (due to the low ionization 
and high column density) at lower energies. Smith et al. (2007) 
analyzed the RGS data and detected two absorption components 
with physical parameters similar (log(£)=2. 14^q ']? 2 and 3.261q'J,7; 
N ff =0.75±g:li and 5.5±^ x 10 21 cm" 2 ) to the ones that we infer, 
strengthening this interpretation. There is also evidence for another, 
higher ionization, mildly relativistic, and variable ionized compo- 
nent in the XMM data. The study of this more extreme component 
is addressed in another paper (Cappi et al., in preparation). 

The observation of highly ionized matter in the core of 
Mrk 509 is in line with its high BH mass and accretion rate. In fact, 
we remind that at the Eddington limit the radiation pressure equals 
the gravitational pull, however the densities of the matter lowers 
with the BH mass (Shakura & Sunyaev 1976). Thus the ionization 
of the material surrounding high accretion rate and BH mass AGNs, 
such as Mrk 509, should be higher than normal. We stress, however, 
that in order to detail the physical parameters of the ionised emit- 
ter/absorber, further long observations are required. 



7 CONCLUSIONS 

The Fe K band of Mrk 509 shows a rich variety of emis- 
sion/absorption components. The XMM-Newton and Suzaku data 
shows evidence for the presence of: 

• a resolved, although not very broad, (a ~0.07 keV) neutral 
Fe Ka line and associated Fe K/3 emission. The width of the line 
suggests that the 6.4 keV line is produced in the outer part of the 
accretion disc (the broad line region or torus emission seem un- 
likely). The measured reflection fraction is consistent in this case 
with the intensity of the line, while a covering factor or column 
density higher than generally observed would be required if the 
line were produced in the BLR or the torus; 

• both the Suzaku and the XMM-Newton data show an excess 
due to ionized Fe K emission. Both datasets show a superior fit 
when a broad ionized line coming from the central parts of the ac- 
cretion disc is considered. The data are inconsistent with narrow 
emission from a distant scattering material at rest, while it can not 
be excluded if the gas is outflowing (v~3500 km s" 1 ) 

• both EPIC-pn and MOS data show an absorption line at ~7.3 
keV, present in the summed spectrum of all XMM-Newton observa- 
tions only. This component confirms the presence of highly ionized, 
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outflowing (v~14000 s -1 ), gas along the line of sight. The com- 
parison between XMM-Newton and Suzaku suggests a variability of 
this component; 

• a hint of an enhancement of variability - both by considering 
the XMM-Newton data alone and by comparison between the two 
data sets - at ~6.7 keV that could be either due to the high variabil- 
ity of the red wing of the broad ionized Fe K line, possibly asso- 
ciated with a variation of the ionisation of the disc, or to a second 
ionized absorption line. 
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